Acoustic wave propagation in an one-dimensional layered system 
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Propagation of acoustic waves in an one-dimensional water duct containing many air filled blocks 
is studied by the transfer matrix formalism. Energy distribution and interface vibration of the air 
blocks are computed. For periodic arrangement band structure and transmission rate are calculated 
analytically, whereas the Lyapunov exponent and its variance are computed numerically for random 
situations. A distinct collective behavior for localized waves is found. The results are also compared 
with optical situations. 
PACS numbers: 43.20., 71.55J, 03.40K 



I. INTRODUCTION 

Propagation of waves in periodic and disordered media 
has been .and, continues to be an interesting subject for 
physicistscTt£l. When propagating in media with inhomo- 
geneities, waves are subject to multiple scattering, which 
leads to many peculiar phenomena such as band struc- 
tures ia periodic media and wave localization in random 
mediaETEi 

The propagation of waves in one dimensional (ID) sys- 
tems has attracted particular interest from scientists be- 
cause in higher dimensions the interaction between waves 
and scatterers is so complicated that the theoretical com- 
putation is rather involved and most solutions require a 
series of approximations which are not always justified, 
making it difficult to relate theoretical predictions to ex- 
perimental observations. Yet wave localization in one 
dimension (ID) poses a more manageable problem which 
can be tackled in an exact manner by the transfer matrix 
methodEllalij. Moreover, results from ID can provide in- 
sight to the problem of wave localization in general and 
are suitable for testing various ideas. Indeed, over the 
past decades considerable progress has been made in un- 
der standjBfi— the localization behavior in ID disordered 
systemsuBll3. However, a number of important issues 
remained untouched. These issues include, for example, 
how waves are localized inside the media and whether 
there is a distinct feature for wave localization which 
would allow to differentiate the localization from residual 
absorption effect without ambiguityE£(Efl. Results from 
the statistical analysis of the scaling behavior in ID ran- 
dom media is not conclusive. Another question could be 
whether the localized state is a phase state which, would 
accommodate a more systematic interpretations. Fur- 
thermore, the study of acoustic propagation in ID ran- 
dom media is relatively scarce. All these motivate us 
to consider wave propagation in ID media further, with 
emphasis on the acoustic wave propagation. 

In this paper we study the problem of acoustic wave 
propagation in an one dimensional water duct contain- 
ing many air blocks either regularly or randomly but on 
average regularly distributed inside the duct. The fre- 
quency band structures and wave transmission are com- 
puted numerically. We show that while our results affirm 
the previous claim that all waves are localized inside an 



ID medium with any amount of disorder, there are, how- 
ever, a few distinctive features in our results. Among 
them, in contrast to optical caseEj, there is no universal 
scaling behavior in the present system. In addition, when 
waves are localized, a collective behaviour of the system 
emerges. We will also show the energy distribution in the 
water duct. 

This paper is organized as follows. The next section 
we explain the the model employed, discuss the transfer 
matrix method and derive relevant formulas. In section 
III numerical results and discussions are given. We then 
summarize the paper in section IV. 



II. MODEL AND METHOD 
A. System setup 

We study the system consisting air blocks inside a wa- 
ter duct. The system is chosen because air filled blocks 
are strong acoustic scatterers. This can be seen as fol- 
lows. The scattering is largely controlled by the acoustic 
impedance which is defined as pc, with p and c being the 
mass density of the medium and the acoustic phase speed 
respectively. The acoustic impedance ratio between wa- 
ter and air is about 3500. This large contrast leads to 
strong scattering, making the system of air blocks in wa- 
ter an ideal candidate for the study of acoustic scattering. 

The ID acoustic system we consider is illustrated by 
Fig. |l|. Assume that N air blocks of thickness Oj (j — 
1, . . . , N) are placed regularly or randomly in a water 
duct with length L measured from the left boundary of 
the duct (LB) . The distance from LB to the left interface 
of the first air block is D. For simplicity while with- 
out compromising generality, in later numerical compu- 
tation we set the thickness of each air block to the same 
value a. The air fraction is clearly j3 = Na/L, the av- 
erage distance between two adjacent air or water blocks 
is (d) = L/N = a//3, and the average thickness of wa- 
ter blocks is (b) = (D + Y^=i b i)/ N - The de g ree of 
randomness for the system is controlled by a parameter 
A in such a way that the thickness of the j-th water 
block is bj = (b)(l + Sj) with Sj being a random number 
within the interval [—A, A]; the regular case corresponds 
to A = 0. An acoustic source placed at LB generates 
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monochromatic waves with an oscillation v(t) = ve~ luJt . 
Transmitted waves propagate through the N air blocks 
and travel to the right infinity. In order to avoid unnec- 
essary confusion, possible effects from surface tension, 
viscosity or any absorption are neglected. 

For convenience, we use the dimensionless quantity 
k(b) to measure the frequency, where k = uj/c is the 
wave number in water blocks. Similarly, k g represents 
the wave number in air blocks. We also define the fol- 
lowing parameters for later use 



(1) 



B. Wave propagation and state vector 

We use the transfer matrix mcthodi for solving the 
wave propagation in the system. Dropping out the 
time factor e~ iuJt , the wave propagation obeys Helmholtz 
equation 



Pm( x ) + k mPm( x ) = °> 



(2) 



in which p m (x) is the pressure field, and the subscript 
to refers to the medium that can be either water or air, 
depending on where x is located. Within any layer (air 
blocks or water blocks), the wave is 



Pm(x) — A m e m + B m e m . 



(3) 



where A r , 



represents the wave transmitted away 



from the source to the right and B m e~ lkmX the wave re- 
flected towards the source. In terms of A m and B m , the 
corresponding velocity field u m , which is another dynam- 
ical variable describing the oscillation of the medium, can 
be calculated as 



u m (x) 



1 



Ibjpr, 
1 



-P'mi x ) 



Jk m x_ g e -ik m xl 



(4) 



where p m refers to the mass density of medium to. 
Define a state vector 



S m (a;) 



S 1 

s 2 




-ik m x 



, (5) 



then p m (x) and u m (x) can be determined by S m (x) 



Pm — S + >S 2 , U Tt 



-(S 1 -^ 2 ). (6) 



We denote the state vector in the j-th air block as Gj (x) 
with two components Gj(x) = Gje tk! > x and G 2 (x) = 
G 2 e~ lkgX and in the j-th water block as ~Wj(x) with 
Wj(x) = Wje lkx and W 2 (x) = W 2 e~ lkx . 

By invoking the condition that the pressure and veloc- 
ity fields are continuous across the interfaces separating 
water and air, one finds 



W^ixj) = JG^xj), Gj(yj) = J- l Wj( yj ), 



with 



J = q 



_i f cosh?] sinhry 
sinh rj cosh r\ 



(7) 



For wave propagation in the j-th air or water layer, we 
have 



Gj(xj) = Ugia^Gjivj), 



(8) 



with 



and 



with 



U g (aj) = 



g ikgdj q 



ik Q (Lj I i 



W j (y j ) = U(b j )W j (x j+1 ), 



(9) 



UQ>i 



ikbj q 



e lkb i 

From these results the transfer matrix Mj for j-th unit 
is 

Mj = JU g {a j )J- 1 U{b j ). (10) 

and the state vectors in water blocks satisfies 

W i _ 1 (a; J ) = M J W i ( a : i+1 ). (11) 

Therefore any two state vectors of the water blocks are 
connected as 

W^-ifoJ = M juj2 W h (x j2+1 ), (12) 

where 

M h , h = M h M jl+1 ■ ■ ■ M h , 1 < ji < j 2 < N. (13) 

From Eq.(|l0|) and (|l^) one can easily prove that all Mj 1} j 2 
are uni- modular matrices, a result of energy conservation. 
We denote 



Mi N EE M = 



Ma M12 



l: ' A - v ' » .VI,; M22 
A simple reduction leads to 

M 2 i=M* 12 , M 2 2=M* 11 . 



(14) 



(15) 



It is clear that the transfer matrix M\,n connects the 
first water block and the last water block. 
Imposing the boundary conditions at LB 

v = -[W^-W 2 }, (16) 



and 

we obtain 
WoCarx) 



pc 



W 2 =0, 



pev 



M n e~ ikD - M* 12 e lkD 



M11 



(17) 



(18) 



The condition in ( p^ ) relates Wo to the LB oscillation 
v, and the second relation ( |l7j ) results from the fact that 
there is no reflection at the rightmost boundary in Fig. 1. 
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C. Parameterization for regular system 



D. Determination of p and u 



In the situation of regular arrangement of air blocks, 
Oj = a, bj = b, dj = d, Mj = M and M = M 1<N = M N . 
According to Eq. (|lo| ) , M acquires the form 



M = 



a P 
[3* a* 



with 



where 



Ikb 



det(M) = 1, 



[cos k g a — zsinfc g acosh2?7], 



(19) 



(20) 



(21) 



(3 = ie sin k g a sinh 2rj, 



are functions of three parameters: k g a, kb and rj. In order 
to express physical quantities in a simple closed form, we 
use the parametrization given inO. We define 

tr(M) = 2cos6», 

then the eigenvalue equation A 2 — tr(Af )A + 1 = yields 
two eigenvalues A = e ±lS . Since tr(M) = a + a* = 2a#, 
and thus an — cos 9, Eq. ( po| ) becomes 

aj - \(3\ 2 = sin 2 9, 

where — Re(a) and aj = Im(a) represent the real 
and imaginary parts of a, respectively. 
Following^, we adapt the parametrization, 



a(6, x) = cos 9 — i sin 9 cosh 2%, 
0(6, X,8)=i e lS sin 9 sinh 2 X , 



(22) 



where 9 and \ m general are complex variables and 5 is 
chosen as kb. Comparing Eq. (|2l|) and ( p2| ) , the following 
relations are established: 

cos 9= coskgd coskb — sinfc g a sinfc6 cosh2?7, 
sin (9 cosh2x= cosk g a sinfcb + sinfc s a coskb cosh 2??, 
sin# sinh2x= sinfc 9 a sinh 27/. (23) 

From these relations one finds that cos 9, sin 9 cosh 2x and 
sin 9 sinh 2% are all real, providing restrictions on the pos- 
sible values of 9 and x- Furthermore, using the definition 
(E2) , it is straightforward to verify that 



[M(d, X ,S)] n =M(n0,x,S) 

a(n9,x) P(n9,x,S) 
P*(n9, x ,6) a*{n6,x) J ' 



(24) 



where n is an integer. 

In this parametrization, matrix M becomes a function 
of 0, Xi an d A direct computation shows that the two 
eigen-vectors of matrix M can be written as 

Xx=f ^ 8hX h ), X 2 =f^t hX V ( 25 ) 
ye 10 sinh x J \ cosn X J 

corresponding to the eigen-values e ±l8 respectively. 



Denoting the pressure and velocity fields on left and 
right boundaries of j-th air block as p^, p^ and Uj, uj*, 
we have 



Pi 



W^x^ + Wf^), 



G](y 3 ) + G 2 ( yj ), 



(26) 



3 

,R 



1 

P c 
1 



PgCg 



G)(y 3 ) - G 2 (y 3 )]. 



The wave vectors Wj-i(xj) and Gj(yj) are determined 
by the following equations: 



W (D) =Mi J _iW i _ 1 (i J ) 

= Mij-iJUg(a j )G j {y j ). 



(27) 



For the case of reg ularly arranged air blocks, by 
Eqs. (|24|), (|26|), and (g7]), p^, pj* and m^, uf are explicitly 
given by 

-a[(N -j + l)6]+ p*[{N- j + l)6]- 



Pf 



pev 



pev 



a(N6) e~ lkD - /3*(N6) e lkD 
a[(N - j)9] e- lkb + 0*[(N - j)6] e lkbl 



i(N9)e- ikD -f3*(N6) e ikD 



(28) 



a[(N-j + l)6]-(3*[{N-j + l)6] 



a(N9) e- lkD ~ f3*(N6) e ikD 
a[(N - 3)6] e- lkb - (3* [{N - j)6] e ikbl 



a(N6) e~ lkD - f3*(N9) e lkD 



E. Band structures 



According to the Bloch theorem, the eigen-modes of 
wave field p and velocity field u in an infinite periodic 
medium can be written as 



p(x) = £(x)e iKx , u(x)=((x)e 



(29) 



where (,(x) and C( x ) are periodic functions satisfying 
£(x + d) = C,(x + d) = ((x) and K is the usual 

Bloch wave number. 
Eq. (§§) implies 



Hid 



W j (x j +d)=W j ^ 1 (a 



M(6, X ,S)W j (x j +d). 



(30) 



Therefore e~ is the eigenvalue of matrix M and equals 
e %s . Substituting Kd = 9 in Eq. (p3|), we obtain the 
dispersion relation 

cos Kd = cos k g a cos kb — cosh 2rj sin k g a sin kb, (31) 
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which describe the band structure. When 

cosk g a coskb — smk g a sinfcfe cosh2r/| < 1, 
the solution for if is a real number. When 

| cosfcgd coskb — smk g a sinfc6 cosh2?7| > 1, 
the solution for K is complex and is of the form 



K=^+i Kl , 
a 



(32) 



where n is an integer in order to satisfy Eq. (plj). The 
frequency ranges within which real solutions for K can be 
deduced define the pass bands, while the ranges rending 
the complex solutions for K determine the stop bands 
(band gaps). 



Rn — 



\/3(N9)\ 



sin 2 NO sinh 2% 



\a(N6)\ 2 l + S in 2 N6 sinh 2 2x' 



and 



rr P V 

o 

pv 

= ~Y 



C + A cos 2N0 + B sin 2N6 
c 



C + Fcos(2N6 
where A, B, C, F, <f> are 

1 



B 



sinh 2x H — cos(2fc£> — 5) sinh 4%, 



sin(2/cL> - S) sinh 2 Xl C = 1 - A, 



(38) 



(39) 



F. Transmission, reflection, energy flow and 
Lyapunov exponent 

Waves propagating through N air blocks will be scat- 
tered many times before they go out. The total transmis- 
sion and reflection rates can be obtained from scattering 
matrix AI S , 



l/t r*/t* 
r/t l/t* 



where t and r are transmission and reflection coeffcients, 
and T = \t\ 2 and R — \r\ 2 define the transmission and 
reflection rates. Based on the previous discussion 



M s = Mi 



N 



and 



T N = 



Rn 



M 



\M 



121 



(33) 



(34) 



Note that in Eq. ( |34| ) the relation Tjv + Rn = 1 is satis- 
fied, which is a consequence of energy conservation since 
the system we considered does not absorb any energy. 

Another consequence of energy conservation is that en- 
ergy flow along the whole duct is a constant. Recall that 
the time averaged energy flow J{x) is calculated as 



J{x) = ^Re[p(x)u*(x)] 



(35) 



In the presence of N air blocks, the flow can be easily 
obtained as 



Jn = 



mi 

2pc 



pv 
~2 



-ikD 



ikD 12 



(36) 



For the case of regularly placed N air blocks we can 
obtain 



T N = 



1 



f 



a(N9)\ 



1 + sin^ N9 sinh^ 2 X 



(37) 



F= VA 2 +B 2 , = tan- 1 (-). 



(40) 



It is well-known that in a random medium waves are 
always localized in space and the localization is charac- 
terized by the Lyapunov exponent (LE) 7. The Lyapunov 
exponent is defined as 



with 



7 = lim (7jv>, 

N— too 



7Ar= 2^ ln( ^ ) - 



(41) 



(42) 



The fluctuation or variance of 7 defined by 

.2 v u._A2^ (43) 



var( 7 ) = lim ((7^} - ( lN ) 2 ) 

N—>oc 



is a quantity which as will be shown below gives impor- 
tant information about the system. 



by 



G. Energy distribution in ID systems 



The time averaged energy density £{x) at x are defined 



£m{x) — 



\'Pn 



2 2 



(44) 



By direct calculations we find that energy density in our 
ID system is piece-wise constant and hence we can sup- 
press the redundant variable x. Energy density in j-th 
water block £f and in j-th air block £ ? are given by 



and 



cw ^ 

~ 2pc 2 



[\w}\ 



\w 



\G 



1 1 2 



\G 



2 1 2 - 



(45) 



(46) 



In the regular case Eq. (45) and (|46j) become 
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ST = 



2pc< 



[M^[(N-j)8]M[(N-j) 



(47) 



and 



si 



where 



H-[M^(N-j + l)6}(J-y 



2 

g c g 



/ = 



xJ" 1 M[(Af - j + 1)6] 



pcv 



11, 



(48) 



Mue- ikD -M* 12 e lkD 
pcv 



a{N6)e- 



ikD 



(3*{N9)e tk 



kD ' 



(49) 



After lengthy but straightforward calculations Eq. ( f47j ) 
and (f48| ) can be further simplified to 



and 



A w + B w cos[2(7V - j) 



£ 9 =C g +A g cos[2(N - 

+B g sin[2(iV-j + l)0] 
= C g +D g cos[2(7V-j + l) 



(50) 



(51) 



where A w , B w , A g , B g , C g , D g and 4> g are defined by 

,,2 



A,. 



pv 



2[C + F cos(2N9 - 0)] 



cosh 2%, 



B„ 



pv 



2[C + F cos(2N9 -4>)\ 



sinh 2\ 



(52) 



(53) 



H. Energy localization and collective behavior 

When waves propagate through media alternated with 
different material compositions, multiple scattering of 
waves is established by an infinite recursive pattern of 
rescattering. Writing p(x) = A(x)e l9 ^ with A(x) = 
\p{x)\ and 9 being the amplitude and phase respec- 
tively, the energy flow J ~ Re[i(p* (x)d x p(x)] becomes 
J - A 2 d x 9. Obviously, the energy flow will come to a 
complete halt and the waves are localized in space when 
A does not equal to zero but phase 9 is constant at least 
by domains. 

From these observations, we propose to use the phase 
behavior of waves to characterize the wave localization. 
Expressing p(x) and u(x) as 



p(x) = A p {x)e l6 ^ x) 



and 



u(x) = A u {x)e l6 ^ x) . 
we construct two unit phase vectors as 
v p = cos6 p e x + sin 8 ' p e y 

and 



(58) 
(59) 
(60) 
(61) 



Physically, these phase vectors represent the oscillation 
behavior of the system. We can plot the phase vectors in 
a two-dimensional plane. 



III. NUMERICAL RESULTS AND DISCUSSION 



A„ 



pv 



Ah[C + F cos(27V0 - <f>)] 



[cosh 2r\ ( 1 — cosh A\) + cos 5 sinh 2r\ sinh 4%] , (54) 



pv 2 sin 5 sinh 2"q sinh 2\ 
9 = 2h[C + Fcos{2N9 - <j))]' 



c = 

9 4h[C + F cos(2N9 - <j>)} 



Aj + B$. 



(55) 



[cosh2?7(l + cosh4x) — cos 5 sinh 2i] sinh 4x] , (56) 



(57) 



The variables C, F, <f> in these equations are the same as 
in Eq. @. 



A. Ordered cases 

In the ordered case, the interference of multiply scat- 
tered waves leads to frequency band structures. For fre- 
quencies located in pass bands, waves propagate through 
the whole system, while for frequencies within a band 
gap, waves are evanescent. 

In Fig. |2|(a), cos Kd versus kb/n from Eq. (|3l]) are dis- 
played for (3 = 10~ 4 (solid curve) and (3 — 3 x 10~ 5 (bro- 
ken curve). The segments of curves between cos Kd = ±1 
(the gray region) give real solutions of Kd and correspond 
to pass bands, whereas those outside of gray region cor- 
respond to forbidden bands. Fig. ||(b) shows the relation 
between the pass/forbidden bands and the fraction of air 
blocks ((3). It is seen that for (3 > 10~ 3 and kb > 3tt, 
the pass bands almost vanish. The band structures for 
various bands are shown in Fig. ||(c) and (d). We see 
that as the volume fraction (3 increases, the width of the 
pass bands decrease and the width of band gaps become 
larger. As expected from the earlier discussion, that air 
blocks are very strong acoustic scatterer in water leads 
to very wide band gaps and narrow pass bands shown in 
Fig. |(c) and (d). 

In Fig. ^, transmission and reflection versus frequency 
for different N are plotted according to Eqs. (p7j) and 
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(pq). Inside a pass band, transmission is significant. The 
transmission rate T/v oscillates as the frequency varies, 
and it has N — 1 peaks as expected from Eq. ([37]). At 
these peaks smN9 = 0, Tn = 1 and Rn — 0. That 
is, at the frequencies waves tunnel through the whole 
system completely, a phenomenon of resonant tunneling. 
Within a band gap, on contrast, although it is not ex- 
actly zero due to the finiteness of the system considered, 
the transmission is significantly prohibited. We have ap- 
proximately Tjy ~ and Rn rj 1. Using Eq. (|39|), energy 
flow Jn as function of the air-block number N for dif- 
ferent frequencies are plotted in Fig. |]. The frequencies 
appearing in diagrams (a), (b), (c), (d) are chosen from 
the first pass band. These diagrams show that the flow 
varies periodically as iV increases. When the frequency 
is low, for instance at kb/ir = 0.02, a simple periodic 
behavior of Jn with long space period is observed. As 
we increase the frequency, the period becomes shorter, 
as illustrated in diagram (b). If we further increase the 
frequency, the oscillation of Jn becomes more rapid and 
the period is smaller than d. However, since TV is a inte- 
ger, the behavior of Jn tends to be more complicated as 
N is enlarged shown in diagrams (c) and (d). 

When frequency enters into the band gap, waves be- 
come evanescent, as indicated in diagrams (e), (f), (g), 
and (h) . Diagram (e) and (h) display the decay behavior 
of Jn at frequencies very close to the gap edges. The be- 
havior of Jn at frequencies near the center of first band 
gap are illustrated in diagrams (f) and (g), showing an 
exponential decay behavior. We see that the decay rates 
in (f) and (g) are much more stronger than in (e) and 
(h), as expected. 

Using Eo^© and @, energy density can be calcu- 
lated. Fig. [| shows the energy density distribution along 
the duct for N — 100 air blocks. On the left panel of 
Fig. H i- e - diagrams (al), (bl), (cl), and (dl), the en- 
ergy density in water blocks are shown, whereas on the 
right panel, i. e. diagrams (a2), (b2), (c2), and (d2), the 
energy density in the air blocks are displayed. In the 
computation, the air fraction is taken as j3 = 10~ 4 . We 
find that for frequencies located in the pass bands, the en- 
ergy density varies periodically along the traveling path. 
For frequencies within the frequency gaps, the energy is 
trapped near the acoustic source, and the energy density 
decays about exponentially along the path. Meanwhile, 
the energy flow is calculated to be close to zero. 

We now look at the oscillation behavior of the blocks. 
As discussed above, zero energy flow leads to a phase 
coherence behavior of the medium. For the purpose we 
consider the behavior of the phase vectors v p and v u de- 
fined in Eqs. (|6(]) and (|6l|). Typical results of phase be- 
havior are shown in Fig.TJ. For further convenience, only 
the phase vectors at the interfaces between air and wa- 
ter are shown. Symbols p L , p R , u L , u R appearing in 
Fig. U denote respectively the phase vectors for the pres- 
sure and the velocity fields on the left and right side of 
the air blocks. We set the vibration phase of LB to 1. 
First we discuss the cases shown by diagrams (a) and (b). 
The frequencies in (a) and (b) are chosen from the first 
pass band. The phase vectors in these two cases point to 
various directions along the duct and v p ■ v u ^ 0, result- 



ing in a nonzero energy flow. The waves are extended 
in the system. In diagrams (c) and (d), the frequencies 
are chosen from the first band gap in which waves are 
evanescent as shown in Fig. B. In this case, all the phase 
vectors of the pressure field are pointing to either n/2 
or — 7r/2, and are perpendicular to the phase vectors of 
the velocity field. The pressure at the two sides of any 
air block is almost in phase. Different from the higher 
dimensional cases in which all phase-.vectors of localized 
fields point to the same directionlllHiS, the present phase 
vectors are constant only by domains. The velocity field 
in neighboring domains oscillate with a phase difference 
7r. At the far end of the sample, however, the phase 
vectors become gradually disoriented, implying that the 
energy can leak out only at the boundary due to the finite 
sample size. We note here that such a phase ordering not 
only exists for the boundaries of the air blocks, but also 
appears inside the whole medium. 

B. Disordered situations 

Unlike in the ordered case for which waves can propa- 
gate through all air blocks if the frequency is located in 
pass band, in a disordered ID system waves are always 
localized. In this section, we study numerically acoustic 
propagation in the disordered case. 

Fig. ^(a) presents the typical results of the transmis- 
sion rate as a function of k(b) for various (3 at a given 
randomness. At frequencies for which the wavelength 
is smaller than the averaged distance between air blocks, 
the transmission is significantly reduced by increasing air 
fraction. 

Fig. 0(b) illustrates the effect of the randomness A 
on transmission for a given air-fraction. For compari- 
son, the transmission in the corresponding regular array 
(A = 0) is also plotted. The gaps are located between 
k(b)/TT = 0.4638 and 1, 1.21 and 2, 2.128 and 3, and so 
on. We find that for frequencies located inside the band 
gaps of the corresponding periodic array, the disorder- 
induced localization effect competes yet reduces the band 
gap effect. To characterize wave localization in this case, 
both the band gap and the disorder effects should be coa- 
sidered, supporting the two parameter scaling theorycJ. 
However, increasing disorder tends to smear out the band 
structures. When exceeding a certain amount, the effect 
from the disorder suppresses the band gap effect com- 
pletely, and there is no distinction between the localiza- 
tion at frequencies within and outside the band gaps. 

Fig. ^ shows that the localization behaviour depends 
crucially on whether the wavelength (A) is greater than 
the average distance between air blocks. When (b)/ A is 
less than 1/4, the localization effect is weak. We also 
observe that with the added disorder, the transmission 
is enhanced in the middle of the gaps. Similar enhance- 
ment due to disorder has also been reported recentlyO. 
Differing fromel, however, the transmission at frequencies 
within the gaps of the corresponding periodic arrays in 
the present system is not always enhanced by disorder. 
Instead the transmission is reduced further by the disor- 
der near the band edges. 
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To further explore the transmission property, we cal- 
culate Lyapunov exponent 7 and its variance var(7) ac- 
cording to Eq. (^|) and ([43]) ■ The sample size is chosen 
in such a way that that it is much larger than the lo- 
calization length and the ensemble average is carried out 
over 2000 random configurations. Fig. presents the re- 
sults for Lyapunov exponent (LE) and its variance as a 
function of k(b)/ir for various randomness. As expected, 
when the randomness is small LE mimics the band struc- 
tures and the variance of LE inside the gaps is small. 
Contrast to the optical caseEJ, there are no double max- 
ima for the variance inside the gap. Rather, the double 
peaks appear in the the allow bands when the system is 
disordered. The double peak feature is more prominent 
in the low frequency bands. When the randomness ex- 
ceeds a certain value, however, the double peaks emerge. 
The higher frequency, the lower is the critical value. For 
example, the double peaks are still visible in the first al- 
low band (c. f. Fig. ||(b)), while there is only one peak 
inside the higher pass bands. Meanwhile, the increasing 
disorder reduces the band gap effect and smears out the 
oscillation in LE, in accordance with Fig. ||. We also 
plot LE versus its variance in Fig. ||. When randomness 
is weak, several branches appear in the LE- variance rela- 
tion. The frequency range of a branch covers a pass band. 
In the optical casa, the frequency range of a branch cov- 
ers a gap insteadEEl. A prominent feature in the present 
system is that when the double peaks in the variance are 
destroyed, the minima of LE correspond to the maxima 
of its variance. This is different from the optical casal3. 

With increasing disorder, we do not observe the gen- 
uine linear dependence between LE and its variance, as 
expected from the single parameter scaling theorycJ, in- 
dicating that the single parameter scaling law may not 
be applicable in the present system. 

Eqs. ( [45| ) and (46) provide the formulas for studying 
energy distribution in any situation. As discussed in the 
regular case, when wave is localized, a kind of ordering of 
the phase vectors appears. From the alignment pattern 
of phase vectors we can know whether waves inside the 
system are localized. For air fraction (3 = 10~ 4 , three 
different cases with different k(b)/ir, A and N are shown 
in Fig. |], [u|[ll]. To isolate the localization effect from the 
band gap effects, we choose the two frequencies located 
in the first allow band: one is in the middle of the band 
and the other is near the lower band edge. 

First, we note that the energy density is constant in 
each individual block. This is a special feature of ID 
classical systems, and can be verified by a deduction 
from Eqs. (45j) and (f46|). From these figures, we ob- 
serve that when the sample size is sufficiently large, waves 
are always localized for any given amount of random- 
ness. When localized, the waves are trapped inside the 
medium, but not necessarily confined at the site of the 
source, unless the band gap effect is dominant. The en- 
ergy distribution does not follow an exponential decay 
along the path- This differs from situations in higher 
dimensionsE3E3. It is also shown that the energy stored 
in the medium can be tremendous. 

These figures also show that for low frequencies and 
when the randomness is small, to trap the waves a large 



number of air blocks is needed. Like in the regular cases, 
when waves are localized, the coherent behavior of the 
medium appears, and the phases are constant by do- 
mains. The phase vector domains are sensitive to the 
arrangement of the air blocks. Moreover, when localiza- 
tion is evident, increasing the sample size by adding more 
air blocks to the far end of the system will not change 
the patterns of the energy distribution and phase vec- 
tors. Therefore the energy localization and the phase 
coherence behavior are not caused by the boundary ef- 
fect. 

Fig. H shows the energy distribution inside the medium. 
Here we choose k(b)/ir — 0.25, which is inside first pass 
band. We see that the energies are localized in both air 
and water blocks. However, unlike in the regular cases, 
the energy does not always decay exponentially along the 
path. In Fig. |l^ we use the same parameters as in Fig. ^ 
except that we choose k(b)/-K = 0.46, which is very close 
to the gap edge 0.4638. We see now that waves can be 
very easily trapped by using only N — 50 air blocks. In 
Fig. [ll] we choose k{b)/ir = 0.25 as used in Fig. ^ but 
increase the randomness to A = 1, waves are trapped by 
N = 150 air blocks, fewer than in the case of Fig. 0. 



IV. CONLUDING REMARKS 

In this paper we studied the propagation of acoustic 
waves in ID layered system consisting of air and wa- 
ter blocks. Both regular and random arrangements of 
air blocks have been studied. We first derived the ba- 
sic formulas for the general situations, and then applied 
these formulas to regular and random cases separately. 
For periodically placed air blocks, the band structures 
were studied and were shown to have large band gaps. 
The transmission, reflection, Luapunov exponent and its 
variance, energy distribution, energy flow and medium 
vibration were also studied. For the case of randomly 
placed air blocks, the results pointed out that waves are 
always confined in a finite spatial region. The disorder 
leads a significant energy storage in the system. It is also 
indicated that the wave localization is related to a col- 
lective behavior of the system in the presence o f mu ltiple 
scattering, also observed for higher dimensionalilllj. The 
appearance of such a collective phenomenon may be re- 
garded as an indication of a kind of classical Goldstone 
modes in the context of the field thcorva. 
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FIGURES 
Figure 1 



□ a 1 a 2 b 2 a 3 b 3 

W | G 1 | W 1 G 2 | W 2 | G 3 | W 3 | ■ ■ 

V =0 x 2 y 2 x 3 y 3 x 4 

FIG. 1. Regular and random arrangements of air blocks 
in the duct. The symbols Gj, Wj denote the air and water 
layers, a-j and bj represent their thicknesses, the j-th air block 
and the water separated by the j-th and j + 1-th air blocks, 
denoted as the j-th water block, is regarded as one unit and is 
labeled as the j-th unit; therefore dj — aj + bj represents the 
thickness of the j-th unit. Xj — D + ^2 j [= ^ di and y.j — Xj+aj 
represent the pisitions of air- water and water-air interfaces. 
For regular arrangement aj — a, bj = b, dj = d = a + b and 
Xj=D + (j - l)d, y 3 = xj + a. 



Figure 2 



DISPERSION RELATION 



BAND STRUCTURE 
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f, = 1e-4 
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FIG. 2. Various properties of periodic system, (a) Dis- 
persion relations for air fractions (3 — 10~ 4 (solid line) and 
/3 = 3 x 10~ 5 (broken line). Segments of curves inside the 
gray region bounded by cos Kd = ±1 correspond to the prop- 
agating waves (allowed bands), whereas those outside of gray 
region represent evanescent wave, (band gaps), (b) Allowed 
bands (gray regions) vs air fraction f3. (c) Band structure for 
(3 = 10~ 4 . (d) Band structure for f3 = 3 x 10~ 5 . Here K is 
the Bloch wave number, d is the lattice spacing. 
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Figure 3 
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FIG. 3. (al)-(a3) Transmission rate versus frequency for 
N = 3, 10, 100. (bl)-(b3) Reflection rate versus frequency 
for N = 3, 10, 100. The volume fraction of air is /? = 10" 4 . 

Figure 4 
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FIG. 4. (a)-(d) Energy flow vs number of air blocks N, 
four frequencies are chosen from first allowed band, (e)-(h) 
Log-Energy flow versus TV in first band gap. The volum frac- 
tion of air is (3 = 10 -4 . 



FIG. 5. Energy density distributions along the duct at four 
different frequencies: kb/n = 0.02, 0.35, 0.4638, 0.9999. Left 
(al)-(dl): Energy density in water blocks. Right (a2)-(d2): 
Energy density in air blocks. Number of air blocks = 100 
and air fraction j3 = 10 -4 . 

Figure 6 
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FIG. 6. Phase vectors at air block boundaries. Total num- 
ber of air blocks is N = 15, and the air fraction is/3 = 10~ 4 . 
Four frequencies are chosen according to Fig. |5_[ that is, 
kb/n — 0.02, and 0.35 in the first allow band and 0.4638, 
0.9999 in the first gap. In the diagrams, the small arrows are 
used to represent the phase vectors. 
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Figure 7 



Figure 9 



(a) Transmission vs Frequency (A=0.3 ) 
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(b) Transmission vs Frequency (p=1e-4 ) 
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FIG. 7. Transmission versus k(b) /n for various air fractions 
at A = 0.3 (a) and different disorders (b). The number of the 
air-blocks is 100. 

Figure 8 



FIG. 9. Energy density distributions along the duct in air 
blocks (a) and in water blocks (b) . Phase vectors at the inter- 
faces for three spatial ranges of the medium are illustrated in 
(c). The unit of energy density is J/m 3 , air fraction (3 = 10~ 4 , 
k(b)/n = 0.25, A = 0.3, N = 1500. 



(a) A=0.008 








' I ' 


: 'i 


' i I 

i 'J 


I S 



0.012 


(d) A=0.008 


0.01 


ft, 


VARIANCE 


f \ 

M 

k 


0.002 







0.04 2 - 5 


(b) 4=0.05 , \ 


2 

0.03 

1.5 

0.02 


/ \ 

f 1 1 1 ' 


°- 0, 0.5 




i U n A 


1 2 3 


0.06 
0.05 




(e) A=0.05 


VARIANCE 


1 

\ 


0.01 





0.05 

j 0.04 
i 

j 0.03 

c 
c 

* 0.02 
0.01 



0.2 

0.15 m 
U 
z 
< 

0.1 E 
> 

0.05 



LYAPUNOV EXPONENT 



LYAPUNOV EXPONENT 



LYAPUNOV EXPONENT 



FIG. 8. Diagrams (a), (b), and (c) show the Lyapunov ex- 
ponent (LE) in broken lines and its variance in solid lines as 
a function of k{b)/7v for three random situations. Diagrams 
(d), (e), and (f) present the plots of the exponent versus its 
variance at the three random cases. Here f3 = 10~ 4 . 
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FIG. 10. Energy density distributions along the duct in air 
blocks (a) and in water blocks (b). Phase vectors at the inter- 
faces for three spatial ranges of the medium are illustrated in 
(c). The unit of energy density is J/m 3 , air fraction f3 = 10~ 4 , 
k{b)/n = 0.46, A = 0.3, N = 50. 
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Figure 11 
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FIG. 11. Energy density distributions along the duct in air 
blocks (a) and in water blocks (b) . Phase vectors at the inter- 
faces for three spatial ranges of the medium are illustrated in 
(c). The unit of energy density is J/m 3 , air fraction f3 = 1CF 4 , 
k{b)/n = 0.25, A = 1,N= 150. 



11 



